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Abstract 

, We study supercurrent conservation for the four-dimensional Wess-Zumino model formulated on 

' the lattice. The formulation is one that has been discussed several times, and uses Ginsparg-Wilson 

^ • fermions of the overlap (Neuberger) variety, together with an auxiliary fermion (plus superpart- 

' ners), such that a lattice version of U{l)ii symmetry is exactly preserved in the limit of vanishing 

I bare mass. We show that the almost naive supercurrent is conserved at one loop. By contrast we 

find that this is not true for Wilson fermions and a canonical scalar action. We provide nonper- 
turbative evidence for the nonconservation of the supercurrent in Monte Carlo simulations. 
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I. INTRODUCTION 



The formulation of supersymmetric field theories on a spacetime lattice is of interest 
because nonperturbative dynamics play an important role in the theory of supersymmetry 
breaking and its transmission to the visible sector of particle physics. Theories such as 
super-QCD and A/" = 4 super- Yang-Mills are of particular interest, but it seems wise to re- 
fine methods using simpler toy models such as the Wess-Zumino model, give the difficulties 
with supersymmetry on the lattice. Hence, we continue our investigations of a lattice for- 
mulation that was studied by several groups a few years ago ■ What we are developing 
here is a methodology for analyzing the extent to which supersymmetry is a feature of the 
low energy effective theory. Since the lattice formulation explicitly breaks supersymmetry, 
this symmetry must be accidental. In fact, it will arise from a fine-tuning of bare lattice 
parameters, corresponding to the ultraviolet definition of the theory. In order to identify 
the supersymmetric point in that parameter space, we must detect the conservation of the 
supercurrent. That is a nontrivial task since the naive supercurrent will mix with other 
operators, due to the explicit violation of the symmetry by the discretization. 

The problems that we face are not by any means unique to the Wess-Zumino model. 
Four-dimensional supersymmetric models on the lattice^ generically require fine-tuning of 
counterterms. This is to be contrasted with lower dimensional theories where lattice sym- 



metries can eliminate the need for such fine-tuning; see lOj for further details. The one 



known four- dimensional exception is pure M = 1 super- Yang- Mills using Ginparg- Wilson 



fermions; the domain wall variety has been the subject of past [13| and recent [9|, Il4j4l8| 
simulations. Clearly we would like to go beyond pure M = 1 super- Yang- Mills, and in 
fact all other models contain scalar fields — which are the source of many difficulties due to 
unwanted renormalizations that cannot be forbidden by symmetries. Recently it was pro- 
posed [iqI that an acceptable amount of fine-tuning could be efficiently performed using a 
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multicanonical Monte Carlo 20|] simulation together with Ferrenberg-Swendsen reweighting 
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in a large class of theories; see also [9|] . In any such program, it is necessary to study 
the divergence of the supercurrent and its renormalization, such as we are doing here for the 
Wess-Zumino model. 



A. Summary of our previous work 

The theory that we study is the four- dimensional Wess-Zumino model, formulated on 



the lattice with overlap (Neuberger) fermions 2J], as well as numerous auxiliary fields. The 



goal of the formulation is to impose the Majorana condition and simultaneously preserve 
the chiral U{1)r symmetry that is present in the continuum in the massless limit. As 



For reviews with extensive references see 
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was shown in our recent work 25|] , preserving this symmetry significantly hmits the number 
of counterterms that must be fine-tuned in order to obtain the supersymmetric continuum 
hmit. In addition to overlap fermions, the lattice formulation has auxiliary fermions (plus 
superpartner fields) that couple to the overlap fermions through the Yukawa coupling, as 



m [26|. It is possible to integrate out the auxiliary fermions (and superpartner fields), and 
when one does this a nonanalytic Dirac operator results for the surviving fermionic field. 
Thus, as has been discussed originally in and at greater length in jj 0| the action is 
singular once auxiliary fields are integrated out. However, as we described in 25|, there is a 
sensible resolution of this singularity by taking the theory to "live" inside a finite box, with 
antiperiodic boundary conditions in the time direction for the fermions. The singularity 
of the Dirac operator that this resolves is related to nonpropagating modes in the infinite 
volume limit; the fact that these are nonpropagating was shown in [^]. However, singularities 

ocalities in the continuum limit, as 



27] . In |25| we measured the degree 



in the Dirac operator raise the spectre of possible non 
was found in gauge theories with the SLAC derivative 
of localization of the Dirac operator following the approach of [28|. We found that while 
there is localization, it is less pronounced than the exponential localization of the overlap 
operator. 

The divergences that need to be cancelled in order to renormalize the lattice theory at 
one-loop turn out to be strictly wave function renormalization. The wave function renor- 
malization of the fermion and the physical scalar match at one loop in the continuum limit 
of the lattice expressions; but, the auxiliary scalar has a mismatched wave function coun- 



terterm. These findings appeared previously in [l| ; thus our work 
those results. 
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was a confirmation of 



B. Plan of this paper 

In Section [TTl we define the Wess-Zumino model both in the continuum and on the lattice. 
We discuss the U{l)ii symmetry of the lattice theory, as well as supersymmetry transforma- 
tions that are a symmetry in the limit of the free theory. Finally we discuss the fine-tuning 
action that must be used to obtain a supersymmetric continuum limit. In Section IIIII we 
describe the supercurrent and an almost naive transcription of it to the lattice. We also 
briefly touch on the form of the renormalized supercurrent in terms of bare lattice opera- 
tors. In Section HVl we perform a one-loop perturbative analysis of the four-divergence of the 
supercurrent. We find that in lattice perturbation theory the supercurrent is conserved. We 
conclude this section by mentioning two-loop diagrams where we expect that the asymmetric 
self-energy of the auxiliary field (at one- loop) will be important, leading to nonconservation 
of the lattice supercurrent. In Section |V] we describe the results of our Monte Carlo simu- 
lations, where we have measured the four-divergence of the supercurrent nonperturbatively. 
We find that the violation of supersymmetry is consistent with contributions beginning at 
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the two-loop order. In Section |VT] we give our concluding remarks. 



II. DEFINITIONS 
A. Continuum 

The Euclidean continuum theory has action 

S = - j d^x \^x^CMx + (i)*U(t) + F*F + F*{m*(t)* + ^0*^) + F(m0 + c/02)|^ 

M = ^ + (m + 2g(j))P+ + (m* + 2/0*)P_. (2.1) 

Our conventions will be {i = 1, 2, 3): 

^± = ^(1±75), C = 7o72. (2.2) 
It can be checked that the action is invariant under the supersymmetry transformations 

5, x = -v^P+(# + F)e- V2P_{^(P* + F*)e, 

6, F = V2e^C^P+x, SeF* = V2e^C^P.x (2.3) 



B. Lattice 



We next discuss the lattice action, which is a special case of the formulations of [l|, [2| ; 



write the lattice action in forms given in 
must be introduced. 



we 



6|. For this, a few lattice derivative operators 



A 
D 



1 - aDw, Dw = ^If^id* + d^) + -ad^d^ 



a 
D 



, + D2 = -{l-A{A^A)-'/') 



(2.4) 
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where 9^ and 9* are the forward and backward difference operators respectively. Note that 
D is the overlap Dirac operator. The lattice action is 



S = —a 



J2 I^^X^CDx + (p*Dl(P + F*F + + F*D2(P* 



--X^CX - - (J^$ + J^*<l>*) 

a a 



+F*{m*(j)* + g*(j)*^) + F{m(f) + 



(2.5) 



As can be seen, the kinetic term for the fermion x involves the overlap Dirac operator D. 
The use of the other operators Di and D2 in the scalar part of the action is a departure from 
what one might do naively, and is the reason for favorable renormalization of the action at 
one-loop. The tilded fields are the linear combinations 



+ $, x = x + x, 



F + J^ 



(2.6) 



The fields J-" and their conjugates are auxiliary fields introduced to allow for a lattice 
realization of the chiral U{1)r symmetry in the m — )■ limit: 

5x = ia'j5 (1 - ^D^ X + iaibX, 5X = iaj^^Dx, 



—3ia(j) + ia 



a 

1 - -D2 
2 



J?* 

2 



(5$ 



-3ia^ + i-a {D24> + F* 



6F = 3iaF + ia 



-D2]F- -Di 
2 J 2 ^ 



2 J,* 



= 3iaF + i-a {D2F + Dl(j)*) 
which takes a particularly simple form on the tilded variables: 

5x = i(^l5Xi ^'P = —2ia(j), 6F = AiaF 
We will only need the supersymmetry transformations of the tilded fields: 

5,0 = ^e^CP+x, 6,^* = V2e^CP^x, 



(2.7) 



(2.8) 



= - V2(P+(Di</) + F)e)is - V2(P_(Di0* + F*)e)fs, 
6,F = V2e^CDiP+x, 6,F* = V2e^CDiP^x 



(2.9) 



This is not a symmetry of the lattice action for g ^ 0, but is a symmetry in the free case. 
Our perturbative analysis in Section |IV] will identify the corresponding conserved current. 
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As noted in [5[, we can integrate out the auxiliary fields X, $, J^, treating the tilded fields 
as constant, to obtain the lattice action:^ 

X ^ 

+F*{m*4>* + g*4>*^) + F{m4) + ^0^) I. (2.10) 



This is the lattice action Eq. (2.14) of with a notation that interchanges Di -ir^ D2, which 
is equivalent to Eq. (2.22) of j2| for the k = case, using the identities^ 

= 75(1-^1)), Tl = l-^D2, D^D = W (2.11) 
The fermion matrix is: 

M = p + mP+ + m*P^ + 2g4)P+ + 2g*4)*P^, ^ = (1 - |d2)"^Di (2.12) 

This way of writing the Dirac operator can be related to the one that appears in ^] by the 
identity: 

{1-^D2)-'D, = {1-^D)-'D (2.13) 
Furthermore we can integrate out the auxiliary fields F, F* to obtain the action 
^ = a'J2[-lx^CMx + l4>*D24>+{m*4>* + 9*4>*'){l-lD2){m^ + g 

This is the action that is used in our numerical simulations. 

When fine-tuning of the lattice action is performed, we must invoke the most general 
lattice action consistent with symmetries. Since we perform our simulations at m 7^ 0, this 
is just the action with all dimension < 4 operators built out of the physical fields, and x- 
We write it here for reference: 

5 = a' 5^ I - Ix'Cip + miP+ + mlP^)x + -J>*D2~^ 

X ^ 

+m2|0|2 + Ai|0|^ + (m202 ^ ^^^3 ^ g^^^*2 ^ ^^^4 ^ ^^^^*3 ^ 
-fC{y,^P+ + yl^*P-)x - x'C{y2~^P^ + yirP+)x\ (2.15) 



^ Integrating out an auxiliary fermion to obtain the fermionic part of this action was previously noted 
m There it was noted that this relates the fermionic action to the one of [2^ by a singular field 
transformation. 

^ We thank A. Feo for explaining this point and providing us with a derivation of these relations. 
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A term linear in cf) has been eliminated by the redefinition )■(/) + c with c a constant. 
The parameters rn^ and Ai are real and all other parameters are complex. Whereas in 
the supersymmetric theory there are four real parameters, in the most general theory we 
have eighteen real parameters to adjust. Holding mi and yi fixed (corresponding to some 
choice of values for and qr in the long-distance effective theory), we have fourteen real 
parameters that must be adjusted to obtain the supersymmetric limit. The counting can 
be alleviated somewhat by imposing CP invariance, so that all parameters can be assumed 
real. Then we have a total of ten parameters. Holding two fixed, we must tune the other 
eight to achieve the supersymmetric limit. Conducting a fine-tuning in an eight-dimensional 
parameter space is a daunting task. 

On the other hand in the limit mi — )■ we can impose the U{1)r symmetry (12. 8p . This 
restricts the action to 

^ = ci'Y.[-\fCIpx + -J>''D^^ + rnl\~^? + M~^\' 

X ^ 

-x''C{y,4>P+ + yl4>'P-)x\ (2.16) 



If we hold yi fixed (corresponding to some choice of qr in the long-distance effective the- 
ory), then only and Ai must be fine-tuned. Conducting a search in a two-dimensional 
parameter space, with both coming from bosonic terms, is manageable. The difficult part is 
that we must extrapolate to the massless fermion limit. Another potential problem is that 
we impose antiperiodic boundary conditions for the fermion in the time direction, but must 
impose periodic boundary conditions for the scalar in order for the action to be single- valued 
on the circle in the time direction. This breaks supersymmetry explicitly by boundary con- 
ditions. At finite mass this should be an effect that can be made arbitrarily small by taking 
the large volume limit. However at vanishing mass, there will be long distance modes that 
will "feel" the breaking due to boundary conditions.'^ Thus it is important that we take 
T ^ 1/ma as m is sent to zero, where T is the number of sites in the time direction. 



What we have seen in 25| is that the one-loop behavior of the theory fl2.10p closely follows 
that of the continuum, so that no new operators are generated at this order. Thus at this 
level of approximation, a fine-tuning of the general lattice action (12.151) is not needed. Due to 
this good one-loop behavior it is of interest to study the original lattice action (12.141) in our 
simulations, without any fine-tuning. By measuring the degree of supersymmetry breaking 
through nonconservation of the supercurrent, we gain information about the higher orders 
and nonperturbative aspects of the lattice theory. 



We thank G. Bergner for raising this point. 
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III. SUPERCURRENT, MIXING AND RENORMALIZATION 



For a general superpotential W{(j)), the supercurrent is 



dW fdW 



(3.1) 

and in our case dW/d(j) = mcf) + gcfp'. There is also a form with the auxiliary field: 

S^ = V2 \^i^P-X + Wl^P^X - F*rP+X - F^P-X] (3.2) 

Because of the supersymmetry breaking on the lattice, this will mix with other operators 
in the same symmetry channel. If the lattice action (Eq. f l2.15p in the massive case and 
Eq. fl2.16p in the massless case) is fine-tuned, then in the long distance effective theory there 
will be a supercurrent that is conserved in the continuum limit. The way to detect the exis- 
tence of supersymmetry in this fine-tuning process is to consider linear combinations of bare 
lattice operators and search for one that has vanishing four- divergence in the supersymmet- 
ric limit, modulo contact terms. We briefiy describe that approach in Section [III Bl below. 
Before doing so, we give an almost naive discretization of the continuum supercurrent (13. ip 
that will turn out to be the conserved supercurrent of the one-loop analysis below. 



A. Almost naive lattice supercurrent 



We formulate this lattice supercurrent in terms of tilded fields: 



dW 

d(f) 



V2 



dW 

ocp 



dW 



liiP-X 



mcj) + gcj) 



(3.3) 



It is almost naive because ^ has been replaced by Di, rather than a more naive prescription 
such as 

1 



(3.4) 



Note that ( 13.30 is the form of the supercurrent with the auxiliary fields eliminated. Thus 
in working with it we should use the Feynman rules corresponding to the action without 
auxiliary fields. In our perturbative calculations below, that will be the action fl2.14p . We 
do not need to use the more general action fl2.15p for the 0{g) calculations that we do, 
since the leading violation of supersymmetry is an 0{g'^) nonsupersymmetric wavefunction 
renormalization for the auxiliary field, as will be discussed below. 
We could also work with a supercurrent containing auxiliary fields: 



D.h^P^X + D,~^*^^P+x - F*rP+X - F^P-X 



(3.5) 
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At finite a, this will lead to slightly different results than (13. 3p . since the lattice equations 
of motion for the auxiliary field are 

-F*= (l-^D2){m4> + g4>^) (3.6) 



In fact we will also consider the form (13. 5 p in our perturbative calculations below. It is a 
convenient choice, since we know crucial facts about the nonsupersymmetric renormalization 
of the auxiliary field from our previous one-loop analysis of counterterms. 

The variation of the action under the lattice version of the supersymmetry transformation 
was given in jsl and is equal to:^ 

6S = ^a^^x^C Lp+(20Di0- ^ ^*p_(20*Di0* - L'l^*^)] e (3.7) 

Note that this is C(a), since 

lim ^ x^CP+{24)Di^ - Di4)'^) = (3.8) 

X 

Also note the presence of g in (13. 7p . In order to see violation of the supersymmetric identity 
dfj,S^ = in the continuum limit, diagrams involving the coupling g must be included. Loop 
diagrams are required, in order to get the divergences that cancel the 0{a) factor coming 
from (13. 7p in the continuum limit. 



B. Renormalized supercurrent 



Having written down the almost naive supercurrent, which will be the subject of all of 
our computations in this paper, we next mention the more general case. That is, the form 
of the renormalized supercurrent, which is expected to be different from the almost naive 
version. At a given engineering dimension, we write down all operators that have the index 
structure of S'^a(x). We denote these as O^^^Jj^ where n/2 is the engineering dimension; thus 
n takes odd values 3, 5, 7, . . ., and the index j labels the different operators of dimension n/2. 
A linear combination of these is the long distance effective supercurrent at lattice spacing a: 

S.ai^)= E (3-9) 

n=3,5,7,... j 

Clearly it will be a demanding task to identify this 5*^ nonperturbatively. However, there is 
no other way to properly fine-tune the action. We must find 5*^ and the point in parameter 
space where it is conserved. Before tacking this problem (in future work), we will (in this 
paper) examine the properties of the almost naive supercurrent. 



^ The variation of the action under a modified supersymmetry transformation was also given in [4|, which 
reduces to the result (|3.7p in the appropriate limits. 
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IV. PERTURBATIVE ANALYSIS 



Here we compute ((9^S'^(a;)(9(0)) for a few different choices of O, with our first choice 
being O = x, which gets an 0{g) contribution at one-loop. As will be seen, this leads 
naturally to a second choice, O = which gets an 0{g^) contribution at one-loop. 

All other choices correspond to higher loop diagrams. As summarized above, we know 
that = ^ Zp at one-loop. These self-energy diagrams can occur on an internal line 
of {d^S^{x)O{0)) , giving for example an 0{g^) contribution when O = (j)*x, at two-loops. 
Because of the mismatch of the Z factors, we can be confident that the self-energies as a 
function of loop momentum are also mismatched, so that the cancellations that occur in the 
continuum theory will not happen. An example diagrams will be presented below. However, 
the point is that we will need to be able to compute two-loop diagrams in order to begin 
fine-tuning the action using the conservation of the supercurrent as a probe. This is beyond 
the scope of the present paper and is left to future work. 

In the perturbative analysis, the following lattice propagators are used: 

-Di{p) + (1 - ^D2{p)){m*P+ + mP_) 



Y,a'e^^-{x{x)x^{0)C) 



iD^ip) + (1 - ID,ip))\ 



m\ 



where 



5^aV^-(0(x)r(O)) = (4.1) 

^ ^^2(P) + (1 - f^2(p))|m|2 



DM = -^^-^E,7.sin(,,a) ^^^^^ 

'[l-2Z,sm'ip,a/2W + Z,sm\p,a) 




l-2E.sin2(p,a/2) 



1 - 2 sin2(p^a/2)]2 + E, sin^p.a) 



(4.3) 



The fermion propagator will be represented by a solid line and the scalar propagator by a 
dashed line. 



A. First choice: O = x 

It is convenient to work instead with the Fourier transform: 

J2e^''-{d'^S,ix)xm (4.4) 

X 

The corresponding diagram at one-loop is Fig. [H Note that this includes the point a; = 
so that we have to worry about contact terms. These would involve the supersymmetric 
variation of x'- 

{Ax{x))6,,o (X {P+{D,4> + F){x) + P-{D,4>* + F*))5,,o (4.5) 
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FIG. 1: The Feynman diagrams for the correlation function that we study. In actuahty, there are 
two such diagrams, one with operator Di^^^P_x on the left-hand side, and one with operator 
07^P+X on the left-hand side. 



At one-loop we found in our previous work {F{x)) = {<p{x)) = 0. So we do not have to 
worry about contact terms at one-loop, due to the absence of tadpoles. 

It is a straightforward calculation to obtain the four correlation functions that we need: 

J2e^'-^d'^{D,h,P-X)a{x)xm) 

X 

^ -g*l{e-'^.- - e^P.«) f ^(D^{q)^,P_S{p - q)P^S(p)CUG{q) (4.6) 
a J [lur 

X 

^ -gUe-^v.a _ ^v.a^ ^ l3_^D,{q)^ ,P^S{p - q) P+ S {p) 0^0 {q) (4.7) 

X 

= -^*l(e-^"'» - e^^'^«) I ^^{^^P^S{p - q)P.S{p)CUG{q) (4.8) 
E^'"''(^M(^V-X)a(^)xMO)) 

X 

^ -^l(e-w« _ eW«) j -0-(7^P_5(p - q)P^S{p)CUG{q) (4.9) 
Notice that these correlation functions have certain common factors. On the left, 

J_(e-^PM« - e^f"") (4.10) 

and on the right 

S{j))C (4.11) 
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a 


Po 


(2^)43(po)i_,oo 


(27r)'^a(po)2-,oo 


0.1 


0.1 


2.0375(2) 


-2.0323(2) 


0.1 


0.2 


4.0700(4) 


-4.0595(4) 


0.01 


0.1 


4.2796(4) 


-4.2789(4) 


0.01 


0.2 


8.5544(9) 


-8.5529(9) 


0.001 


0.1 


6.5529(7) 


-6.5513(7) 


0.001 


0.2 


13.1009(13) 


-13.0979(13) 



TABLE I: Values of the integrals for m = 1, at various values ol pq and a. 



Thus when we sum them, we can factor out these bits and the remainder has to be the thing 
that cancels. We now check whether the four-divergence of the almost naive supercurrent 
vanishes. We have computed the above integrals numerically for various values of a, having 
in mind the a — )■ extrapolation. 

If the naive supercurrent is to work, then with p = {po, 0, 0, 0) 

^iPo),.,.p = J -0^,{D,iqhoP-Sip-q)P-)aMQ) 

^iPo)2-,ap = I -0^,iloP+S{p-q)P^UGiq) (4.12) 

the quantity &{pq)_ has to vanish in the continuum limit. It is not hard to show that a 
naive continuum limit of this expression does vanish, provided 

= ^2, feS^/'^ = V2ma (4.13) 

i.e., for the proper coefficients of the almost naive supercurrent. If &{po)_ „^ vanishes in the 
one- loop continuum limit, it would also follow that 

&{Po)^,^p = j ^|&r^*Pi(g)7oP+^(p-g)P+)«/3G(g) 

+6f/')*a-^(7oP-^(p-g)P+).^G(g)| (4.14) 

would vanish as well. 

When one looks at the specific values of the integrals, (Table H]) what one finds is that we 
should set the coefficients to the naive values of fl4.13p . (Of course the overall normalization 
of \/2 is a matter of convention.) That is, the almost naive supercurrent is conserved at one 
loop. 

This result suggests that at one-loop there is an argument that the quantity we are 
computing is determined by the free theory and the supercurrent is not renormalized at 
one-loop. 
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FIG. 2: The simpler correlation function that must vanish in the continuum limit. 



B. Second choice: O = (t)*x 

That the forgoing claim is true can be seen from the fact that when we amputate the 
propagator S{p) from the diagram Fig. [1] (i.e., the fermion propagator that is not in the 
loop), we obtain the diagram of Fig. O However, the latter diagram is just the one that we 
would obtain in the free theory from evaluating (5^S'^(x)(0*x)(O)). Since in the free theory 
the variation of the action under supersymmetry is zero, there should be a conserved current. 
What we have just found is that the almost naive lattice supercurrent is that supercurrent 
to within the accuracy of our numerical evaluation of the integrals. 



C. A more naive discretization 



The behavior we have just observed is to be contrasted with what happens in the free 
theory if we use Wilson fermions and naive scalars so that the action is 

S = -Y^ a\^X^C{Dw + m)x + (p*d^d;^ + F*F + m{F(j) + (4.15) 

X 

where Dw is defined in Eq. (12. 4p and we have specialized to a real mass m. Also, we will 
use the symmetric difference Dirac operator (13. 4p in the discretization of the supercurrent 
(EH), ,^^^5, yielding 



Ds{p) = -- 7m sin(p^a) (4.16) 



in momentum space. Corresponding to (I4.15p . after integrating out the auxiliary field, there 
will be fermion propagator 

^^^^ ^ sin^(p,a) + (m + f sin\p,a/2W ^ ' 
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a 


Po 


(27r)4j(po)3_,oo 


(27r)4a(po)4-,oo 


0.1 


0.1 


7.772923 


-2.076717 


0.1 


0.2 


15.540440 


-4.148735 


0.01 


0.1 


65.802(7) 


-4.4650(4) 


0.01 


0.2 


131.599(13) 


-8.9252(9) 


0.001 


0.1 


623.85116 


-6.77184 


0.001 


0.2 


1247.69800 


-13.53882 



TABLE II: Values of the more naive integrals for m = 1, at various values of pq and a. 
and scalar propagator 



It is easy to check that 
transformation, such as 

5a = -V2P+{Ds<P + F)e - V2P-{Ds<P* + F*)e, 

6,F = V2e^CDsP+x, S,F* = V2e^CDsP-X (4.19) 

in spite of the fact that this is a free theory. I.e., it does not have the behavior (13. 7p seen in 
the discretization that is the main focus of this paper. 
The integrals that are the counterparts of (14.12^ are: 

f d'^Q 

^(Po)4-,„/. = I -0^,boP+Sw{p-q)P-UGM (4.20) 

It can be seen from Table [III that setting (I4.13P will not work: there is a significant lat- 
tice artifact and the naive supercurrent is not conserved. This shows the benefits of the 
formulation which is the main focus of this article. 

D. Higher dimensional operators 

Another interesting thing we have looked at is 

/ ^A^^(^)^^i^)^oP-S{p - q)P-)^rMq) (4.21) 

which corresponds to one of the dimension 9/2 operators that would appear in (13. 9p . The 
power of a could possibly be overcome by a 1/a coming from the extra Di{q). In fact what 



^E^sin2(p^a/2) +m2 ^^'^^"^ 
the action (I4.15P is not invariant under a lattice supersymmetry 
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FIG. 3: A diagram where the nonsupersymmetric renormahzation of the auxihary field will play 
a role. The dashed/dotted line is the mixed {F(f)) propagator. The dotted line is the (FF*) 
propagator. 

we find is that the integral is linearly divergent. This indicates that higher dimensional 
operators could play an important role in the renormalized supercurrent when we go to 
higher orders. Obviously we would prefer to avoid this possibility, since there is an infinite 
number of such operators. However, we do not at this stage see a reason for excluding them 
from the sum (13. 9p . 

E. Beyond one-loop 

The result of the one-loop investigation is that there is a conserved supercurrent in the 
continuum limit, when we work to 0{g). This agrees with the fact that the renormahzation 
at one loop only yields a nonuniform wavefunction renormahzation for the auxiliary field: 
Zp ^ = (in the continuum limit). This cannot affect correlation functions at 0{g) 
since the wavefunction renormahzation of the auxiliary field would have to appear on the 
internal line of a diagram, which necessarily implies a factor of g"^ to occur. For example, 
the diagram in Fig. |3] would be sensitive to the value of Zp, and is O^g"^). Note that for 
this analysis we have switched to the supercurrent involving the auxiliary field, Eq. (13. 2p . 
In the nonperturbative analysis that we discuss next, we find further evidence that the 
nonconservation of the almost naive supercurrent begins at Oi^g"^). 

V. NONPERTURBATIVE ANALYSIS 

Here it is convenient to take the spatial transform so that we instead work with 

E ^) = E ^tS,o.{t, x) = dfQ^it). (5.1) 

X X 
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ma 


9 


TJi^a^ Ai lattice 


0.1 


0.1 


— — 8^ X 32 


0.1 


0.1 


0.01 0.01 8^ X 32 



TABLE III: Parameters of the Monte Carlo simulations that we have performed. 



We evaluate 

Cit) = J2{d',S,{t,^)Om (5.2) 

X 

Since S^a is fermionic, an odd number of x fields must appear in nonvanishing correlation 
functions with dtQaif)- The case that we will consider is C = 'X^C. Because of cluster 
decomposition, C{t) will fall off with t exponentially, governed by the mass meg of the 
lightest state created by O. As we tune the action, meg and the fall-off with t will change. 
Thus we could mistake a decrease in C{t) for an improvement of supersymmetry when it is 
really an increase in meg. Similarly, we could think we have worsened supersymmetry when 
in fact all we did was to decrease meg. Clearly we need a way to normalize C [t) in order to 
cancel off this exp(— t rries) behavior. For this reason we look instead at the ratio 

_ IEx(^M^..(^.x)(x^C),(0))| 

IEx(x.(t,x)(Fc^),(o))| ^'-'^ 

and will set a = /3 = 0. 

Here we use Monte Carlo simulations to nonperturbatively measure the ratio (15.31) using 
;he almost naive supercurrent (13. 3p . The simulation method is rational hybrid Monte Carlo 



29|, and the runs were performed on Compute Unified Device Architecture (CUDA) enabled 
graphics processing units (Nvidia GeForce GTX 285, GTX 480 and Tesla C1060), using code 
that we developed and tested in our previous work. We have measured the autocorrelation 
time to be approximately 12 molecular dynamics time units for ma = 0.1 bare fermion mass, 
with a coupling g = 0.1. The length of the simulation was 5,000 molecular dynamics time 
units and we sample at each 5 time units. Errors in the ratio function R{t) are computed by 
jackknife analysis with data blocked into 5 samples each. The runs are summarized in Table 
mil As can be seen we also consider the case of the fine-tuning action (I2.16P except that we 
allow for a bare fermion mass m and set yi = g. Results for R{t) in each case are displayed 
in Figs, m and [5l It can be seen that at large times, where the long distance theory should 
be obtained, R{t) ^ 0{g'^). This is consistent with nonconservation of the almost naive 
supercurrent beginning at two loops or higher. It can also be seen that the action (I2.14p 
gives a significantly smaller value for R{t) than the action (I2.16P — again, a fermion mass 
term has been added to the latter. This shows how the formulation (I2.14p has a particularly 
small violation of supersymmetry, even at higher orders. 
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FIG. 4: The ratio ([O]) for ma = 0.1, 5 = 0.1. 
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FIG. 5: The ratio 1^ for ma = 0.1, g = 0.1, m^a^ = 0.01, Ai = 0.01. 



VI. CONCLUSIONS 



We found that at one loop (9^5'^ (x) 0(0)) = 0. This was true without any tuning of the 
lattice action, and provided the almost naive supercurrent is used. The numerical results 
were explained by the fact that the one-loop diagram is related to a free theory diagram, 
and so must vanish. We showed that this result does not hold if a more naive discretization 
is used. Next we discussed two loop diagrams where we do not expect the cancellations to 
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hold, since they are sensitive to the mismatch in self-energies that was already found in our 
previous study of one-loop counterterms. We look forward to presenting numerical results 
for the two loop diagrams in a forthcoming paper. Finally, we provided nonperturbative 
results with the almost naive supercurrent. It was seen that the nonconservation of the 
supercurrent is consistent with contributions beginning at two loops. Another direction for 
future research is the fine-tuning of the action together with the search for the renormalized 
supercurrent, which will have the more general form fl3.9p . Investigations in this direction 
are in progress. 
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